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ABSTRACT 


Rigid-body  slamming  was  experimentally  investigated  at  the 
Naval  Ship  Research  and  Development  Center  by  dropping  one  flat 
circular  aluminum  plate  model  and  five  cone-shaped  aluminum 
models  with  t.,.11  deadrise  angles  (up  to  15  deg)  from  various 
elevated  positions  above  a  calm  water  surface.  This  report 
presents  the  test  results  and  compares  them  witii  theory  and 
to  those  for  two-dimensional  wedge-shaped  models. 

ADMINISTRATIVE  INFORMATION 

This  investigation  was  carried  out  as  part  of  a  general  stud}’  on 
slamming.  The  work  was  funded  jointly  by  Subproject  S46-06X,  Task  1707 
(Hydrofoil  Hull)  and  by  JSESi’O.  However,  the  publication  of  this  report 
was  funded  by  JSESPO  alone. 

INTRODUCTION 

The  two-dimensional  impact  test  of  a  rigid  flat-bottom  model  indi¬ 
cated  that  the  maximum  impact  pressure  is  nowhere  near  the  theoretical  in¬ 
finitely  large  hydrodynamic  pressure  nor  near  the  theoretical  acoustic 
pressure.1  The  cushioning  effect  of  the  compressible  air  trapped  between 
the  impact  body  and  the  water  surface  reduces  the  maximum  impact  pressure 
to  about  one-tentli  of  the  acoustic  pressure.  However,  the  nature  of  the 
trapped  air  phenomenon  is  not  very  stable.  Much  more  air  is  trapped  for 
the  impact  of  a  flat  bottom  and  a  1-deg  wedge  than  for  wedges  with  deadrise 
angles  of  4  deg  or  higher. 

Tests  of  elastic  models  verified  the  fact  that  the  pressure 
generated  by  the  impact  is  affected  by  the  vibratory  movement  of  the  im¬ 
pact  surface  and  that  it  can  be  separated  into  rigid  body  impact  pressure 
and  interacting  pressure.1  This  dynamic  interaction  is  closely  related  to 
the  hydrodynamic  phenomenon  rather  than  to  the  acoustic  phenomenon. 

Past  studies  have  therefore  demonstrated  that  in  the  two-dimensional 
impact  of  rigid  and  elastic  bodies  (1)  water  can  be  treat  ?d  as  an 


^References  are  listed  on  page  40. 
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incompressible  fluid  regardless  of  the  size  of  the  deadrise  angle  of  the 
wedge,  (2)  trapped  air  must  be  taken  into  consideration  for  small  deadrise 
angles,  and  (3)  the  structural  response  to  impact  can  be  treated  as  the 
impact  of  a  deformable  b<  y  on  an  incompressible  fluid,  with  or  without 
trapped  air. 

As  a  con  equence  of  these  findings,  questions  naturally  arose 
regarding  the  three-dimensional  effect  of  slamming  and  led  to  the  series 
v ,  drop  tests  of  cone-shaped  models  discussed  herein.  In  order  to  compare 
two-  ar  i  three-dimensional  impact  phenomena,  drop  tests  of  a  two-dimensional 
flat-bo.com  model  were  repeated,  and  test  data  of  wedge-shaped  models  wore 
reanalyzed.  A  theoretical  investigation  on  slamming  of  cone-shaped  bodies 
was  also  performed;  results  3re  preset  •■i  sewhe.-e.  “ 

The  objective  of  the  present  study  is  therefore  to  investigate  the 
three-dimensional  effects  of  slamming  from  the  drop  tests  of  wedges  and 
cones.  To  fulfill  this  objective,  the  trapped  .  ■ r  phenomenon  and  pressure 
time  histories  are  compared  for  these  ..cries  of  ti. ts. 

The  measurement  of  transient  pressure  requires  instrumentation  of  a 
specialized  nature.  Appendix  A  of  the  report  is  devoted  to  a  description 
of  instrumentation,  features  of  circuit  designs,  calibration  techniques, 
and  operational  procedures.  To  utilize  the  General  Electrr:  Computer  Time¬ 
sharing  System,  a  small  program  was  written  to  fit  the  test  results  by  the 
method  of  least  squares.  This  program  is  presented  in  Appendix 

Drop  tests  of  cone-shaped  models  were  performed  at  the  Center  during 
Fiscal  Year  1969. 


METHOD 


DESCRIPTION  OF  ^’-TS  FOR  STUDYING 
THREE-D1MENS10NAI  EFFECTS 

The  six  16- i n . -diameter  aluminum  models  consisted  of  a  flat  circular 
plate  with  O-deg  deadrise  angle  and  five  cone  shapes  with  deadrise  angles 
of  1,  3,  6,  10,  and  15  deg.  The  impact  surfaces  of  the  models  were 
machine-finished  to  ensure  smoothness  and  true  conical  configuration. 
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Tests  were  conducted  in  a  25-  by  IS- ft  rectangular  tank  which  was 
filled  with  water  to  a  depth  of  7  ft  2  in. 

Details  of  the  test  assembly  are  shown  in  figure  1.  The  releasing 
mechanism  consisted  of  a  solenoid  attached  to  a  cross  beam  by  an  adjustable 
steel  rod  which  can  be  raised  and  lowered  for  the  proper  drop  height  of 
the  model.  The  solenoid  was  equipped  with  a  hook  for  hanging  the  model. 
When  the  solenoid  was  activated,  the  hook  was  very  quickly  released  and 
the  model  fell  freely  in  the  vertical  direction.  The  total  drop  weight 
of  the  assembly  shown  in  the  figure  was  05  lb. 

The  drop  heights  (deJV,  .  ike  distance  between  the  vertex  of  the 
cone  and  the  water  surface)  ranged  >  .o  15  in.  in  5-in.  increments. 

Pressures,  accelerations,  vertical  displacement,  and  velocities  of  the 
me.  ng  model  were  record-  '.  In  addition,  10-mm  high-speed  underwater 
movies  were  taken  to  ?  i>  the  trapped  air  phenomenon  during  and  after  the 
impact.  The  film  spovd  varied  up  to  about  1000  frames  per  sec  for  the 
study . 

The  instrumentation  system  used  for  the  experiment  consisted 
essentially  of  quartz-crystal  transducers,  charge  amplifiers,  d-c  ampli¬ 
fiers,  and  a  tape  recorder.  The  validity  of  the  pressure  measurements  of 
the  complete  recording  system  was  tested  electronically  and  mechanically; 
the  system  was  also  calibrated  by  means  of  an  underwater  explosion.  The 
results  indicated  that  the  entire  recording  system  had  the  ability  to  pick 
up  and  record  any  high- frequency  acoustic  pressure  existing  during  the 
impact  of  the  falling  body  with  the  water  surface."5  The  detailed 
description  of  instrumentation  is  given  in  Appendix  A. 

ADDITIONAL  DATA  IJTIbIZLi) 

Reanalyzed  test  data  on  rectangular  flat-bottom  and  wedge-shaped 
models  were  utilized  for  comparison  with  results  of  the  cone-shaped 
models.  Details  of  procedure  and  analysis  are  available  in  Reference  4. 


difference  Bin  wi-:i:n  test  SETUPS  OF  CONE- 
SHAPED  MODELS  AND  WEDGE-SHAPED  MOPE  1. 5 


The  rectangular  flat-bottom  ami  wedge-shaped  models  were  tested  in 
the  25  by  15  ft  rectangular  tank  with  two  parallel  walls  to  ensure  two- 
dimensional  flow  conditions.  These  two  parallel  walls,  made  of  steel 
plates,  were  constructed  to  span  the  length  of  the  tank  and  to  extend 
from  18  in.  above  the  water  surface  to  the  full  tank  depth.  They  were 
rigidly  connected  to  the  tank  floor  and  sides  and  were  separated  by  a 
distance  equal  to  the  model  length  of  20  1/2  in.  plus  a  small  amount  of 
clearance.  These  parallel  walls  had  open  ends  to  permit,  free  flow  of  the 
surface  wave  around  the  tank  during  the  drop  test. 

The  cone-shaped  and  circular-plate  models  were  tested  in  the  same 
tank  but  without  these  parallel  walls  so  that  the  flow  of  air  and  water 
was  not  restricted. 

THREE-DIMENSIONAL  EFFECT  ON  TRAPPED  AIR 

High-speed,  10-mm  underwater  movies  of  the  drops  of  wedge-shaped 
and  cone-shaped  models  revealed  that  only  the  0-  and  the  1-deg  deadr.ise- 
angle  configurations  of  both  shapes  trapped  a  considerable  amount  of  air 
and  pushed  some  of  it  into  the  surface  layer  of  water.  Most  of  the  air 
was  not  trapped  at  the  instant  of  impact  by  either  shape  with  deadrise 
angles  of  3  deg  and  higher.  During  the  impact,  the  higher  the  deadrise 
angle  of  both  the  wedge  and  the  cone,  the  clearer  and  cleaner  was  the  view 
of  the  impact  surface.  However,  for  the  same  deadrise  angle,  the  wedge 
trapped  more  air  than  did  the  cone.  This  phenomenon  is  illustrated  by  the 
underwater  photographs  shown  in  Figures  2  and  3. 

The  fact  that  the  three-dimensional  (cone-shaped)  body  allowed  the 
trapped  air  to  escape  earlier  than  did  the  two-dimensional  (wedge-shaped) 
body  was  evident  in  comparisons  of  the  detection  of  air  trapped  between 
the  impact  surface  of  the  falling  model  and  the  water  surface;  see 
Figures  4  and  5.  The  method  used  to  detect  trapped  air  is  given  in 
Reference  4.  The  record  for  the  wedge  impact  showed  that  only  after  the 
first  positive  pulse  of  the  impact  pressure  was  completely  over  did  the 
trapped  air  appear  partly  to  have  escaped  and  partly  to  iiave  been  pushed 
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into  the  surface  layer  of  water.  However.,  the  record  for  the  cone  impact 
showed  that  the  trapped  air  began  to  escape  or  to  be  pushed  through  the 
water  surface  during  the  first  positive  pulse  of  the  impact  pressure.  Be¬ 
cause  of  this  difference,  their  impact  pressure  time  histories  were 
affected,  especially  for  the  impact  of  the  1-deg  wedge  and  cone.  This  will 
be  discussed  further  in  the  next  section. 

THREE- DIMENSIONAL  EFFECT  ON  IMPACT  PRESSURE 


The  pressure  time  histories  of  the  two-dimensional  and  three- 
dimensional  impacts  were  affected  because  of  differences  in  the  time 
required  for  trapped  air  to  escape.  These  differences  can  be  observed  in 
Figures  4  to  7.  The  shape  of  the  pressure  time  histories  were  very  much 
alike  except  for  the  1-deg  models,  where  the  three-dimensional  effect  on 
the  escape  time  of  the  trapped  air  reduced  the  cushioning  effect  of  the 
air.  Thus,  the  appearance  of  the  impact  pressure  time  histories  of  the 
1-deg  cone  is  similar  to  those  of  cones  with  higher  deadrise  angles,  but 
the  appearance  of  the  impact  pressure  time  histories  of  the  1-deg  wedge 
closely  resembles  those  of  the  two-dimensional  flat-bottom  model. 

The  maximum  impact  pressures  for  both  shapes  of  models  are  plotted 
on  the  log-log  charts  shown  in  Figures  8  and  9.  As  shown  in  these  figures, 
the  data  recorded  for  the  wedge-shaped  models  were  much  less  repeatable 
than  those  of  the  cone-shaped  models  because  of  improvements  in  the  drop 
mechanism  and  the  instrumentation  system  used  for  the  cone-shaped  models. 
The  pressure-velocity  relations  were  obtained  empirically  from  test  data 
by  fitting  straight  lines  on  the  charts  with  a  general  form  of 

p  =  k  7 

where  k  and  n  are  the  arbitrary  constants  and  n  has  limits  of  1  <  n  <  2. 

The  method  just  described  provides  equations  to  estimate  the  maximum 
impact  pressures  of  wedges  and  cones  penetrating  a  water  surface;  they  are 
summarized  as  follows: 


lm  Two-Dimensional  Wedges  and  Flat  Bottom: 


a.  Flat  bottom  (20  in.  wide): 


At  center: 

I’max  =  0-443  p  V2 

X 

6  11/16  in,  off  center: 

Pmax  ~  0.376  p  V“ 

! 

b.  1-deg  wedge: 

| 

Away  from  keel: 

pmax  =  °-5i5  p  V2 

I 

j 

! 

c.  3-deg  wedge: 

i 

Away  from  keel: 

pmax  =  O'387  P  V2 

i 

d.  6-deg  wedge: 

< 

a 

Away  from  keel ; 

Pmax  =  °-387  P  V2 

e.  10- deg  wedge: 

Away  from  keel: 

Pmax  =  0.186  p  V2 

i 

f.  15-deg  wedge: 

i 

Away  from  keel: 

Pmax  =  °-10S  P  V2 

Three-Dimensional  Cones  and  Circular 

Plate: 

i 

a.  Circular  plate: 

At  center; 

pmax  -  °-34S  P  V1-96 

1 

4  in.  off  center: 

Pmax  =  °'304  p  V1-90 

i 

b.  1-deg  cone: 

Away  from  vertex: 

Pmax  =  2-387  P  V1’60 

e.  3-deg  cone: 

Away  from  vertex: 

Pmax  =  °'778  P  V1'83 

d.  6-deg  cone: 

Away  from  vertex: 

Pmax  =  °’366  P  V1-84 

! 

1 

e.  10- deg  cone: 

Away  from  vertex; 

pmax  =  °-160  P  V1-91 

;  i 

f.  15-deg  cone: 

!  j 

Away  from  vertex: 

Pmax  =  °’129  P  V1’70 

i 

I  i 

1  i 
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The  formulas  for  the  wedges  arc  from  Reference  1  with  slight 
revisions.  The  formulas  for  the  flat  bottom  were  obtained  from  a  recent 
test,  which  is  believed  to  be  better  than  the  previous  tests,  by  using  the 
method  of  least  squares  given  in  Appendix  B.  The  formulas  for  the  cones 
and  circular  plate  were  also  obtained  from  a  recent  test  by  the  method  of 
least  squares. 

If  test  results  of  the  cones  and  circular  plate  are  forced  to  fit 
into  velocity-square  relations  as  were  those  for  two-dimensional  wedges  and 
flat  bottom  models,  the  formulas  for  them  are: 


a.  Circular  plate: 

At  center: 

4  in.  off  center: 

b.  1-deg  cone: 

Away  from  vertex: 

c.  3-deg  cone: 

Away  from  vertex: 

d.  6-deg  cone: 

Away  from  vertex: 

e.  10-deg  cone: 

Away  from  vertex: 


p 

1  max 
p 

rmax 

p 

1  max 

p 

1  max 

p 

*  max 

pmax 


0.320  p  V2 
0.258  p  V2 

1.160  p  V2 

0.562  o  V2 

0.273  p  V2 

0.134  p  V2 


f.  15-deg  cone: 

Away  from  vertex: 


p  =  0.072  p  V 
1  max 


Comparison  of  test  results  for  the  wedges  and  cones  indicated  that 
the  maximum  impact  pressures  were  generally  lower  for  the  three-dimensional 
impact;  see  Table  1.  However,  because  of  the  reduced  cushioning  effect  of 
trapped  air  (see  and  compare  pressure  time  histories  of  1-deg  wedges  and  ( 

1-deg  cones  shown  in  Figures  6  and  7) ,  the  maximum  impact  pressures  of  the 
1-deg  cone  were  higher  than  those  of  the  1-deg  wedge.  For  comparison, 
the  equations  just  listed  and  theories  of  two  wedges  and  cones  ’  are  also 
plotted  in  Figures  8  and  9.  Since  the  cushioning  effect  of  the  trapped  air 
is  gradually  reduced  as  the  deadrise  angle  is  increased,  the  theoretical 
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and  the  experimental  maximum  impact  pressures  at  the  same  impact  velocity- 
are  much  closer  in  magnitude  for  the  higher  deadrise  •  angles  of  wedges  and 
cones  than  for  the  lower  deadrise  angles. 

TABLi:  1 


Comparisons  of  k  Values  in  Wedge  and  Cone  Impact 
Pressure  empirical  Formulas 


Deadrise 

Angle 

deg 

k  values 

2 

in  formula  p  =  kpV 

'  max 

Cone 

Wedge 

Cone/Wedge 

Flat  (center) 

0.320 

0.443 

0.72 

1 

1.160 

0.  515 

2.25 

3 

0.562 

0.887 

0.64 

6 

0.273 

0.387 

0.71 

10 

0. 134 

0. 186 

0.72 

15 

0.072 

0.  103 

0.70 

SUMMARY  OF  FINDINGS 


Previous  studies  on  slamming  at  the  NSRDC  Department  of  Structural 
Mechanics  involved  tests  for  two-dimensional  cases.  The  present  study  ex¬ 
tended  the  investigation  to  the  three-dimensional  aspect  and  studied  its 
effects  on  slamming. 

Drop  tests  were  conducted  on  six  aluminum  models  having  conical 
impact  surfaces  with  respective  deadrise  angles  of  0  (circular  plate),  1, 
3,  0,  10,  and  IS  deg  (see  Figure  1).  On  the  basis  of  this  series  of  ex¬ 
perimental  investigations,  the  following  conclusions  are  drawn: 

1.  During  the  rigid-body  impact  of  the  circular  plate,  the  first 
positive  pulse  of  the  impact  pressure  occurs  at  the  time  when  the  trapped 
air  begins  to  escape  (Figure  5).  Only  the  circular  plate  and  the  1-deg 
cone  trap  considerable  amounts  of  air;  cones  with  deadrise  angles  of  3  deg 
and  higher  do  not  trap  very  much  air  (Figure  3). 

2.  The  maximum  impact  pressures  are  generally  about  30  percent 
lower  for  the  cone  than  for  the  wedge  (Figure  8) . 
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3.  Because  the  trapped  air  escaped  more  quickly  after  impact  by  t lie 
cone  than  by  the  wedge  shape,  the  pressure  time  histories  for  the  1-dcg 
cone  closely  resemble  hydrodynamic  impact  rather  than  impact  with  cushion¬ 
ing  air  (see  Figure  7).  In  contrast,  the  pressure  time  histories  for  1-deg 
wedge  resemble  impact  with  cushioning  air  (see  Figure  G) .  Thus,  the  maxi¬ 
mum  impact  pressures  of  the  1-deg  cone  are  higher  than  those  of  the  1-deg 
wedge  (Figure  8) . 
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Rectangular  Plate 


6-Deg  Wedge 


1-Deg  Wedge 


10-Deg  Wedge 


3 -Deg  Wedge 


15-Deg  Wedge 


Figure  2  -  Detection  of  Trapped  Air  during  Drop  Tests  of 
Rigid  Wedge-Shaped  Models 
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Circular  Plate 


6-Deg  Cone 


3-Dug  Cone.  1 3 -Dug  Cone 

Figure  3  -  Detection  of  Trapped  Air  during  Drop  Tests  of 
Rigid  Cone-Shaped  Models 
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Figure  5  -  Sample  Record  Showing  Time  Lag  between  Pulse  of 
Maximum  Impact  Pressure  and  F.scape  of  Trapped  Air  for 
Impact  of  Three-Dimensional  Rigid  Circular  Plate 
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figure  7  -  Sample  Records  Taken  during  0-lncli  Drop  Tests 
of  Rigid  Cone-Shaped  Models 
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Figure  7  (Continued) 
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Figure  7  (Continued) 


Run  180:  15-Deg  Cone 
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MAXIMUM  IMPACT  PRESSURE  W  PS  I 


MAXIMUM  IMPACT  PRESSURE  IN  PS  I  y  MAXIMUM  IMPACT  PRESSURE  IN  PS  I 


Figure  9  -  Experimental  Results  of  Maximum  Impact 
Pressure  Caused  by  Water  Impact  of  Rigid  Cone- 
Shaped  Models 
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MAXIHUM  IMPACT  PRESSURE  IN  PS  I 


APPENDIX  A 
INSTRUMENTATION 


GENERAL  DESCRIPTION 

The  general  layout  of  the  instrumentation  array  is  shown  in  the 
functional  block  diagram  of  Figure  Al. 

The  body  force  measuring  channels  consisted  of  three  piezoelectric 
(PE)  pressure  gages,  one  linear  PE  accelerometer,  and  four  identical  charge 
amplifiers.  The  charge  amplifier  output  volt,  ge  signals  were  amplified  and 
recorded  on  Channels  2  through  S  of  tie  magnetic  tape  recorder. 

The  vertical  displacement,  Channel  6,  consisted  of  a  linear  po- 
tentiometrical  displacement  gage  and  a  displacement  box  which  included  a 
constant-current  gage  drive  and  an  amplifier  to  drive  the  tape  recorder. 
This  box  also  included  a  differentiator  circuit  for  obtaining  a  vertical 
velocity  signal  which  was  recorded  on  Channel  7.  A  secondary  calibration 
system  referenced  to  the  basic  physical  quantities  of  distance  (in  feet) , 
time  (in  seconds)  and  velocity  (in  feet  per  second)  was  included  in  this 
same  unit  (labeled  "displacement  box"  in  Figure  Al) . 

A  water  contact  sensor,  Channel  8,  and  a  1-kc  sinusoidal  timing 
reference,  Channel  1,  completed  the  eight-channel  datum  input  to  the 
Ampex  CP- 100  magnetic  tape  recorder. 

A  calibration  switch  box  allowed  a  sequential  electrostatic  charge 
calibration  of  the  PE  body  force  measuring  channels  before  the  model  was 
released.  This  unit  also  activated  the  drop  solenoid  which  released  the 
model  at  the  end  of  the  sequential  charge  calibration. 

The  data  were  recorded  at  60  IPS  and  played  back  at  the  end  of  the 
run  at  1  7/8  ips  onto  a  7-in.  CEC  oscillograph  for  immediate  on-the-spot 
inspection  and  analysis. 

GENERAL  FEATURES  OF  CIRCUIT  DESIGNS 

Displacement  and  Velocity  Measurement  System 

The  displacement  and  velocity  measurement  system  is  illustrated  in 
the  systemic  block  diagram  of  Figure  A2.  The  general  features  of  the 
circuit  designs  of  each  systemic  section  are  described.  The  theoretical 
aspects  will  be  discussed  in  a  later  section. 
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Figure  A1  -  General  Layout  of  Instrumentation  Array 
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Figure  A2  -  Displacement  and  Velocity  Measurement  System 


Displacement  Potentiometrical  Transducer:  The  vertical  displacement 

* 

measurement  transducer  consisted  of  a  1 OK- ohm,  10-turn  potentiometer.  The 
body  of  this  transducer  was  affixed  to  the  stationary  frame  of  the  "drop" 
apparatus.  Mounted  on  the  shaft  potentiometer  was  a  sheave  pulley  with 
serrations  to  receive  a  light  beaded  chain  that  connected  to  the  model. 

The  beaded  chain  passed  up  over  this  pulley  and  was  tensioned  by  means  of 
a  small  lead  weight. 

As  shown  in  Figure  A3,  two  field  effect  transistors  (PUT)  were 
connected  in  series  with  a  2.5K-ohm  resistor  which  shunted  the  displacement 
potentiometer.  The  lower  FHT  served  as  a  constant-current  (cc)  generator 
to  energize  the  transducer.  The  lK-ohm  rheostat  in  this  cc  circuit  was 
adjusted  so  that  when  the  beaded  chain  was  moved  exactly  1  ft  (thus  rotating 
the  shaft  of  the  potentiometer) ,  the  resultant  voltage  change  on  the 
slider  was  precisely  0.5v.  The  upper  F13T  circuit  was  adjusted  so  that  with 
the  model  just  touching  the  water  surface,  the  voltage  on  the  slider  was 
precisely  zero.  The  operational  follower  amplifier  provided  the  low  im¬ 
pedance  required  to  drive  the  differentiator  circuit  in  the  displacement 
velocity  convertor. 

Calibration  Signal  Generator:  The  circuitry  for  the  calibration 
signal  generator  is  shown  schematically  in  Figure  A4.  A  high-gain 
operational  comparator  (741C  on  left  side)  was  coupled  by  a  resistance 
train  to  an  operational  integrator  (741C  center).  The  comparator  abrupt ly 
changed  state  when  its  two  input  terminals  (Terminals  4  and  S)  were  equal. 

A  square  wave  with  an  amplitude  of  14  VP  was  integrated  to  provide  a 
symmetrical  triangle  at  Terminal  10  of  the  integrator  (741C  center).  The 
amplitude  of  this  triangular  wave  was  about  6  VP.  This  amplitude  was 
trimmed  to  precisely  2  VP  by  the  lOK-ohm  potentiometer  shown.  A  unit-gain 
follower  (741C  right  side)  provided  low- impedance  isolation.  The  timing 
period  was  precisely  1  sec,  ±  sec.  Trimming  adjustment  is  provided 


K  =  10^  ohm. 
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by  the  5G0k-ohm  po. „.itiometer  in  the  resistance  train.  The  zero  line  and 
amplitude  stability  of  this  total  circuit  were  excellent,  being  in  the 
order  of  a  part  in  5000  or  better. 

Output  Card  Tab  13  led  to  Input  Card  Tab  4  on  the  parabolic  wave 
generator  shown  in  Figure  AS. 

Parabolic  Wave  Generator:  The  parabolic  wave  generator  is  shown 

schematically  in  Figure  A5.  This  unit  is  best  described  as  a  highly 

sophisticated  operational  integrator.  The  integrator  proper  is  shown  on 

the  left  side  of  the  illustration.  The  right-hand  circuitry  was  a  high- 

accuracy  zero-line  stabilizer  which  constituted  the  principal  difference 

between  this  and  conventional  operational  integrating  circuits.  The  two 

* 

RC  networks  (3.5M/.047  p f,  330K/.47  pf)  acted  as  a  negative  feedback  ab¬ 
sorption  filter,  which  removed  the  dynamic  output  from  the  integrator. 

The  zero- line  or  d-c  component  passed  without  attenuation  through  the 
right  hand  741C  amplifier  which  had  a  gain  of  times  10.  The  "clean" 
output  of  this  amplifier  was  fed  to  positive  input  Terminal  5  of  the  inte¬ 
grator.  This  became  100  percent  "negative"  feedback  because  of  the  phase 
reversal  in  the  X10  amplifier.  The  resultant  stability  factor  for  the 
zero-line  was  about  one  part  in  50,000. 

Displacement  Velocity  Converter:  The  displacement  velocity  con¬ 
verter  consisted  of  a  conventional  operational  differentiator  shown 
schematically  at  the  upper  left  side  of  Figure  A6.  The  displacement  signal 
entered  on  Card  Tab  2.  The  sensitivity  of  the  converter  had  a  fine  ad¬ 
justment  control,  the  IM-ohm  potentiometer  shown.  This  control  was  ad¬ 
justed  as  follows.  The  parabolic  wave  train  with  an  amplitude  of  0.5  VP 
was  imposed  on  Card  Tab  2.  This  represented  a  displacement  of  1  ft.  At 
Terminal  10  of  the  741C  (upper  left)  differentiator,  the  issue  had  to  be  a 
triangular  wave  train  with  a  period  of  precisely  1  sec  and  an  amplitude 
of  4  VP.  This  represent'd  a  velocity  of  8  fps.  In  other  words,  the 
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106  ohms  and  pf  =  10~6 


farad. 
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Figure  A5  -  Parabolic  Displacement  Wave  Generator 
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Figure  A 6  -  Displacement  Velocity  Converter 


28 


triangular  output  wave  from  the  differentiator  had  to  be  twice  the  exact 
replica  of  the  2-VP  triangular  wave  that  was  impressed  on  the  parabolic 
wave  generator.  The  factor  of  two  results  from  the  practical  requirement 
of  a  1/2  sec  time  constant  in  the  parabolic  wave  generator.  This  repli¬ 
cation  was  the  heart  of  the  accuracy  of  the  entire  velocity  system.  After 
this  fine  adjustment  was  completed,  the  true  displacement  signal  from  the 
pcitentiometrical  transducer  (which  had  been  calibrated  to  yield  0.5  VP/ft) 
could  be  switched  through  this  same  differentiator  circuitry  and  thus 
read  out  true  velocity. 

In  order  to  fit  the  amplitude  scale  of  the  tape  recorder,  the  4-VP 
triangular  wave  was  passed  through  the  1/8-attenuator,  thence  through  the 
x2  operational  amplifier  on  the  right  side.  Now  1  VP  corresponded  to  a 
velocity  of  8  fps.  The  lower  x2  operational  amplifier  served  to  convert 
the  1/2  VP  displacement  signal  to  1  VP  or  lv  per  foot  of  model  trav  1. 

Water  Contact  Sensor 

The  water  contact  sensor  was  designed  to  indicate  the  exact  instant 
that  the  model  touched  the  water  surface.  As  shown  in  Figure  A1 ,  a  30K-ohm 
resistor  was  connected  to  the  5-kc  calibration  oscillator.  Two  small 
foil  plates  were  attached  on  the  bottom  of  the  model;  these  were  centrally 
located  but  insulated.  One  plate  connected  to  the  30K-ohm  resistor  and 
also  went  to  an  amplifier  for  recording  on  "mag"  tape.  The  other  plate 
connected  to  the  signal  ground.  When  the  water  was  contacted,  these  two 
plates  were  thereby  essentially  short -circuited  and  the  amplitude  of  the 
5  kc  wave  train  was  reduced  to  approximately  one- tenth.  This  provided  a 
simple  but  effective  indication  of  the  time  of  contact  with  the  water. 

CALIBRATION  TECHNIQUES 
Body  Force  Calibration 

The  body  force  measuring  transducers  consisted  of  three  PE  pressure 
gages  and  one  linear  PE  accelerometer.  Four  identical  charge  amplifiers 
were  utilized  to  electronically  convert  the  charge  signals  to  voltage 
signals  for  amplification  and  recording. 
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Because  a  PE  transducer  inherently  issues  a  charge  (q)  signal 
rather  than  a  voltage  (v)  signal,  a  basic  charge  calibration  rather  than  a 
voltage  calibration  was  used  to  calibrate  the  entire  system.  This  cali¬ 
bration  was  based  on  the  basic  relationship  (Coulomb's  Law): 

Q  =  CV 

The  charge  calibration  signal  was  generated  by  using  the  S-kc  oscillator 

shown  in  Figure  A1  and  driving  this  voltage  into  the  calibration  input  of 

the  respective  charge  amplifiers.  The  calibration  inputs  consisted  of 
*• 

100-pf  capacitors  between  the  calibration  and  signal  inputs  of  the 
amplifier.  Thus  a  charge  calibration  signal  (q)  was  generated  which  was 
the  combination  of  the  S-kc  oscillator  (voltage  v)  and  the  100-pf 
capacitor  (capacitance  c) . 

The  amplitude,  and  thus  the  peak  value,  of  this  charge  calibration 

signal  was  adjusted  precisely  by  the  microvolter  shown.  The  sensitivities 

of  the  individual  PE  gages  were  given  either  in  10~12  coulombs  per  pound- 

per-squarc-inch  or  per  g.  Now  a  known  calibration  of  the  approximate  unit 

of  measure  peak-to-peak  could  be  impressed  on  each  respective  channel.  An 

*  * 

illustration  of  this  conversion  of  units  is  presented  below, 
charge  calibration  signal  (pcb  P-P) 

- — - — — - =  calibration  value  (psi  P-P) 

gage  sensitivity  (pcb/psi) 

where 

gage  sensitivity  =  X  pcb/psi 
C  =  100  pf 
NX 

v  =  - v  rms 

2^2 

1  v  rms  =  2V2  V  (P-P) 


*  -12 
pf  =  10  farad. 

**  -12 

pcb  P-P  =  10  coulombs  peak  to  peak. 
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Theoretical  Aspects  of  the  Free- 
Fall  Calibration 

The  free-fall  calibration  system  was  of  the  secondary  type.  It 
automatically  produced  analog  voltages  that  were  referenced  to  the  basic 
physical  units  of  distance  (in  feet)  time  (in  seconds),  and  velocity  (in 
feet  per  second) .  The  method  permitted  on- record  comparison  between  the 
actual  motion  of  a  dropped  body  and  the  electronically  simulated  motion 
of  a  body  in  free  fall  under  ideal  conditions. 

A  functional  wave  generator  (TIMER  1  in  Figure  A2)  delivered  a 
triangular  wave  train  which  was  accurately  symmetrical  about  its  zero  line 
and  which  had  a  periodicity  of  precisely  1  sec.  This  wave  train  had  a 
single  peak  amplitude  of  2v,  or  2  VP. 

The  triangular  wave  train  is  impressed  on  an  operational  integrator 
(Item  2).  The  issue  from  the  integrator  is  a  symmetrical  parabolic  wave 
train  also  with  a  periodicity  of  precisely  1  sec  and  having  a  single  peak 
amplitude  of  0.5v  or  0.5  VP. 

With  vhe  CAL-OPERATE  switch  in  the  CAL  position,  the  parabolic 
wave  train  was  .impressed  on  an  operational  differentiator  (Item  4  (d/dt) 
on  Figure  A2) .  The  issue  from  the  differentiator  was  a  replica  of  the 
original  triangular  wave  train  with  an  amplitude  of  4  VP.  After  passing 
through  a  fixed  attenuator  (x  1/8)  and  ar.  amplifier  (x2),  the  triangular 
wave  train  had  an  amplitude  of  1  VP  which  represented  and  corresponded  to 
a  linear  velocity  wave  train  with  a  single  peak  amplitude  of  8  fps . 

Simultaneously,  the  parabolic  wave  train  passed  through  another 
amplifier  (x2)  and  issued  with  an  amplitude  of  1  VP  which  represented 
and  corresponded  to  a  free-fall  displacement  wave  of  1  ft. 

The  significance  and  meaning  of  this  parabolic  wave  train  requires 
some  discussion.  The  physical  situation  being  electronically  simulated 
was  that  of  an  imaginary  acceleration  field  with  an  acceleration  force  of 
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1  g  continuously  directed  toward  a  line  or  plane.  Thus  if  an  object  is 
held  1  ft  from  this  plane  and  then  dropped,  it  will  fall  under  a  constant 
accelerating  force  g  until  it  passes  through  the  plane.  At  this  instant, 
the  force  encountered  is  exactly  reversed,  whereupon  the  velocity  de¬ 
creases  to  zero,  the  object  "falls"  back  through  the  plane,  and  in  1  sec 
again  reaches  the  position  from  which  it  was  initially  dropped.  In  the 
assumed  absence  of  frictional  or  any  other  forces,  the  object  follows  this 
parabolic  displacement  pattern  indefinitely. 

The  accelerating  force  acting  on  the  object  just  described  is  shown 
in  Figure  A7a  as  a  symmetrical  square  wave.  Figure  A7b  depicts  the 
velocity  pattern  which,  of  course,  is  a  symmetrical  triangular  wave. 

Figure  A7c.  presents  the  displacement  pattern  which  is  the  parabolic  wave 
train  under  discussion. 

A  single  parabolic  arch  of  this  displacement  wave  train  is  shown  in 
Figure  A8.  The  concept  represented  by  this  illustration  is  as  follows: 
an  object  with  a  velocity  of  8  fps  directed  upward  along  OH  emerges  at  the 
origin  0.  This  object  continues  to  travel  upward  at  decreasing  velocity 
until  the  height.  11  is  reached,  at  which  point  its  velocity  is  zero. 
Immediately  thereafter,  the  body  falls  from  H  through  the  origin  0  at  a 
velocity  which  increases  linearly  with  time. 

The  parabolic  arch  or  curve  is  a  true  time-history  of  the  displace¬ 
ment  of  the  body  from  the  time  of  its  emergence  at  the  original  0  until  it 
falls  through  the  same  point  on  exit.  This  curve  serves  to  spread  the 
motion  on  a  convenient  time  scale  to  show  the  detail  of  the  solely 
vertical  excursion.  There  is  no  horizontal  motion,  and  none  should  be  in¬ 
ferred. 

The  classical  formulation  for  a  parabola  of  this  type  and  kind  is: 

2 

x  =  2  py 

where  p  is  the  sizing  parameter  (i.e.,  P/2  is  the  distance  from  the  focus 
to  the  vertex)  and  2p  spans  the  width  of  the  curve  in  a  sector  cut  through 
the  focus. 

Since  the  displacement  s  of  a  falling  body  is  described  by  the  ex¬ 
pression: 
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S  =  1/2  gt“ 

the  classical  formulation  is  modified  by  adjusting  the  nomenclature  and  by 
moving  the  origin  0  to  the  left-hand  side  as  shown  in  the  illustration. 
This  modification  and  translation  is  accomplished  as  follows: 

Step  1.  Let  x  =  t,  y  =  h,  and  p  =  -1/g.  Then 

t2  =  (-2/ g)h 

Step  2.  Move  origin  0  to  (-T,  -II).  Then 
(t  -T)2  =  (-2/g)  (h-H) 

where  h  is  the  instantaneous  height  above  the  origin  0, 

H  is  the  maximal  height  of  h, 
t  is  the  instantaneous  time  notation,  and 

T  is  the  particular  value  of  t  at  midpoint  under  the  vertex. 

Step  3.  Set  limits  0  <  t  <  2T.  Then 

<*  -n2  u»  t .  ».2T  ■  <-2«  <M» 

arid  for  this  particular  calibration  system:  H  =  1  ft,  g  =  32  ft/sec,2  and 
T  =  1/4  sec.  Then 

(t  -1/4)  lim  t  _  0+l/2  =  "  f6  Oi-l) 

OPERATIONAL  PROCEDURES 

The  system  was  connected  as  shown  in  the  general  layout  of  the 
instrumentation  array,  Figure  A1 .  Calibration  and  gain  setting  values  for 
the  various  model  shapes  and  drop  height  conditions  were  preset.  These 
values  were  readjusted  or  updated  as  experience  warranted  prior  to  each 
run  or  drop.  The  potentiometrical  transducer  was  adjusted  so  that  at  the 
water  level  it  read  precisely  zero  on  the  DVM  provided  for  this  purpose. 
The  model  was  then  elevated  to  the  desired  drop  height  and  suspended  by 
the  drop  solenoid.  The  Ampex  CP-100  magnetic  tape  recorder  was  set  up  to 
record  at  60  ips . 
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With  these  initial  setup  conditions  met  and  the  proper  warmup  con¬ 
ditions  for  each  systemic  unit  satisfied,  the  following  steps  have  to  be 
performed  in  sequence: 

1.  Turn  the  displacement  box  switch  to  the  CAL  position. 

2.  Turn  the  tape  recorder  to  RECORD . 

3.  After  approximately  5  sec,  turn  the  displacement  box  back  to  the 
OPERATE  position. 

4.  Next,  slowly  rotate  the  drop  switch  on  the  calibration  switch 
box  from  Position  1  to  Position  7  -  this  sequentially  calibrates  the  PE 
gage  channels  and  then  releases  the  model. 

5.  Approximately  5  sec  after  the  model  drops,  stop  the  tape 
recorder. 

6.  Rotate  the  drop  switch  on  the  calibration  switch  box  back  to 
Position  1. 

7.  Turn  the  displacement  box  switch  back  to  the  CAL  position. 

8.  Rewind  the  tape  and  playback  at  1  7/8  in/sec  on  the  7  in.  CEC 
osci 1 lograph . 

9.  Analyze  the  record  and  adjust  gain  settings  and  CAL  values  as 
required  for  the  next  run. 
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APPENDIX  B 


A  FORTRAN  IV  PROGRAM  TO  FIT  TEST  DATA  BY 
THE  METHOD  OF  LEAST  SQUARES 


Assume  that  the  test  data  points  can  be  fitted  with  the  equation 

,  ,,n 
p  =  k  V 


which  may  be  rewritten  as 

log  p  =  log  k  +  n  log  V 


This  equation  is  a  form  of 


Y  =  a  +  a.  X 
o  1 


If  there  are  N  numbers  of  test  data  points,  then  by  the  method  of  least 
squares, 

N  N 

Z  VN  ■  %  N  *  “l  Z 
1  1 

N  N  N 

Z  XN  VN  Z  X»  *  al  Z  XN 

l  l  l 


These  two  equations  give 


a  = 
o 


IN  IN  IN  IN 

Z  VN  Z  V  Z  XN Z  XN  VN 


r  n 


NZ  <-Z 


L  1 


Thus 


k  =  log  * 

If  the  test  data  points  are  forced  to  fit  into  the  equation 

p  =  k  V2 

For  N  numbers  of  test  data  points, 

N  N 

£  \  £  <vvn2> 

i  i 

k  = -  =  - 

N  N 

The  computer  program  for  this  simple  mathematical  operation  is  as 
follows: 

10  DIMENSION  X(.100),Y(100),P(100),V(100),Z(100) 

12  READ,NCASES 
15  D0  50  NC=1,NCASES 
20  READ.N 
30  D0  10  1-1, N 
40  10  READ, PCI) ,V(I) 

50  XN=N 

60  D0  20  1=1, N 

70  Y ( I ) =AL0G  (PCI)) 
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